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Abstract 

It has been recently shown that the observed baryon number may originate 
at the electroweak phase transition, provided that the Higgs boson and the 
lightest stop are sufficiently light. In this work, we perform a detailed anal- 
ysis, including all dominant two-loop finite temperature corrections to the 
Higgs effective potential, as well as the non-trivial effects proceeding from 
the mixing in the stop sector, to define the region of parameter space for 
which electroweak baryogenesis can happen. The limits on the stop and 
Higgs masses are obtained by taking into account the experimental bounds 
on these quantities, as well as those coming from the requirement of avoiding 
dangerous color breaking minima. We find for the Higgs mass nih < 105 GeV, 
while the stop mass may be close to the present experimental bound and 
must be smaller than, or of order of, the top quark mass. These results pro- 
vide a very strong motivation for further non-perturbative analysis of the 
electroweak phase transition, as well as for the search for Higgs and stop 
particles at the LEP and Tevatron colliders. 
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1 Introduction 



Recently, there has been renewed interest in the idea of generating the observed baryon 
number |IJ at the electroweak phase transition A physical realization of this 

idea seems to be possible only by going beyond the Standard Model formulation. In 
fact, if no fundamental standard model singlet scalar field is present in the low energy 
effective theory, electroweak baryogenesis seems to be naturally induced only if the Higgs 
mass is sufficiently light, and if new light scalar particles, with strong couplings to the 
Higgs field are present in the theory |J. Both effects tend to induce a strongly first 
order electroweak phase transition [01-0, which allows to preserve the generated baryon 
number. Moreover, new CP-violating phases, beyond those present in the Standard 



Model, seem to be necessary to fuel the baryon number generation [10 



Minimal low energy supersymmetric extensions of the Standard Model fulfill the 
three above mentioned requirements. A light CP-even Higgs naturally appears in these 
models. Indeed, an upper bound, of order 125 GeV can generally be set on the lightest 



CP-even Higgs mass [ 12 ]-[ 14 1 . The superpartners of the top quark, whose coupling to 
the Higgs is of order of the top quark Yukawa coupling, provide the additional scalars 
which help in inducing a relatively strong first order phase transition ||15|- ||17|| . For 
this to happen, the lightest stop should have a mass of order of, or smaller than, the 
top quark mass |I5[-|P5[. Moreover, new CP- violating phases, induced by the super- 



symmetry breaking parameters associated with the left-right stop mixing, are naturally 



present in these theories ||11||. A very recent analysis shows that an explanation of the 



observed baryon number of the universe requires CP-violating phases of order one and 



chargino and neutralino masses not much larger than 200 GeV P^fl . Similar results were 
obtained |25j by different analysis (see also Refs. |26], |27|, p8|). 

The above results present a strong case for the search of Higgs and top-squark 



particles at LEP and the Tevatron colliders ||29|| . It is therefore of the highest interest to 



define, in the most accurate possible way, the exact region of parameter space consistent 
with electroweak baryogenesis in the rrih-mj plane. As has been first shown in Refs. [^, 



23fl , two loop effects are very important to define the correct boundary of this region. In 
Ref. [|3U] it was suggested that, in the case of zero mixing in the stop sector, the phase 
transition can be sufficiently strong for Higgs masses as high as 100 GeV, if the stop 
mass is inside a narrow band, of order 150 GeV. 

In this work, we define the allowed parameter space, including all dominant two- 
loop corrections to the Higgs effective potential and allowing a non-trivial mixing in the 
stop sector. We show that stop mixing effects highly increase the region of stop masses 
consistent with electroweak baryogenesis. Since light stops may induce color breaking 



minima, which should be taken into account in the analysis |18| , the finite temperature 
effective potential along the squark and Higgs directions is computed. We identify the 
regions of parameter space in which the physical vacuum is stable within the thermal 
evolution of the universe, as well as those in which a two-step phase transition, via a 
color breaking minimum, can take place. Finally, we also define the region of parameter 
space in which the system is driven to the physical vacuum state at high temperatures, 
but this state becomes unstable at low temperatures. We define the boundaries which 
separate the regions consistent with the different above mentioned possibilities and show 
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that, even in the case of absolute stability of the physical vacuum, stops as light as 100 - 
120 GeV may be consistent with electroweak baryogenesis, for Higgs masses in the range 
80-100 GeV. 

In section 2 we summarize the most important properties of the effective potential 
at finite temperature. The expressions of the effective potentials for the stop and Higgs 
fields, with all dominant finite temperature two-loop corrections included, are given in 
appendix A and B, respectively. In section 3 we define the different regions of parameter 
space to be studied in this article. In secion 4 we present our results. We reserve section 
5 for our conclusions. 

2 Finite temperature effective potential 

In order to preserve the generated baryon number after the electroweak phase transition, 
a large vacuum expectation value of the Higgs field at the critical temperature is required 



where v = ((f)) and v(0) ~ 246 GeV. The phase transition must be hence strongly first 
order. The implications of the above bound may only be obtained by a detailed knowl- 
edge of the electroweak phase transition. In the Standard Model, it has been shown 
that the two-loop improved effective potential leads to values of the order parameter 
and the critical temperature in good agreement with those obtained through a lattice 
non-perturbative analysis |7|, || |3i~f . In this case the requirement of Eq. ( |2.1|) demands a 
value of the Higgs mass which has been already ruled out by experimental searches. In 
this work, we shall present a two-loop perturbative analysis for the case of the Minimal 
Supersymmetric Standard Model (MSSM), in the expectation of a complete numerical 
simulation of this model. 

Following the methods developed in Ref. M, we compute the four dimensional finite 
temperature effective potential keeping all dominant two-loop finite temperature cor- 
rections induced by the Standard Model gauge bosons as well as by the light squarks 
and Higgs particles. As has been first shown in Ref. p2|| , the two-loop corrections to 
the Higgs effective potential in the Minimal Supersymmetric Model are dominated by 
loops involving the light stops and the SU(3) gauge bosons, and strongly depend on 
the QCD coupling. 

Although we reserve the technical details of the effective potential calculation for 
appendix A, we summarize here the most important properties of the two-loop finite 
temperature effective potential for the Higgs field. First of all, at one loop, the dominant 
corrections to the effective potential come from the gauge bosons, the top quark and 
its supersymmetric partners. We shall work in the limit in which the left handed 
stop is heavy, itlq > 500 GeV. In this limit, the supersymmetric corrections to the 
precision electroweak parameter Ap become small and hence, this allows a good fit to 
the electroweak precision data coming from LEP and SLD []. Lower values of rag make 

1 A lighter left-handed stop may also be consistent with precision electroweak measurements p^ ]. 
However, this requires sizeable values of the stop mixing mass parameter, which tend to suppress the 



i 





2 



the phase transition stronger and we are hence taking a conservative assumption from 
the point of view of defining the region consistent with electroweak baryogenesis. The 
left handed stop decouples at finite temperature, but, at zero temperature, it sets the 
scale of the Higgs masses as a function of tan (3. For right-handed stop masses below, 
or of order of, the top quark mass, and for large values of the CP-odd Higgs mass, 
m,A 3> Mz, the one- loop improved Higgs effective potential admits a high temperature 
expansion, 



m 2 (T) , 



E SM 3 + (2i\Q 



K+n ra (T) 

12tt 



3/2 



A(T) 



= + . . . 
(2.2) 



where 



[ 2 ^ 2 ' 1 ^h + sm 2 (3(l-A 2 t /m 2 Q 



9 



6 



/ (l_±; t ™23;:...'7 ' (2.3) 



is the finite temperature self-energy contribution to the right-handed squarks (see ap- 
pendix A), g s is the strong gauge coupling, N c = 3 is the number of colours and 
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is the cubic term coefficient in the Standard Model case. 

Within our approximation, the lightest stop mass is approximately given by 



m\ ~m 2 v + 0.15M| cos 2(3 
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where A t = A t — \ij tan (3 is the stop mixing parameter and m t = h t <p sin P/y/2. As was 
observed in Ref. the phase transition strength is maximized for values of the soft 
breaking parameter mfj = —Hj(T), for which the coefficient of the cubic term in the 
effective potential, 



E ~ ^SM + 



h 3 t sin 3 (3 (l- A 2 /m 2 c 
4V2vr 



3/2 



which governs the strength of the phase transition, 



v(T c ) 



AE 



(2.6) 



(2.7) 



can be one order of magnitude larger than E$m PI. In principle, the above would 



allow a sufficiently strong first order phase transition for Higgs masses as large as 100 
GeV. However, it was also noticed that such large negative values of m\j may induce the 



presence of color breaking minima at zero or finite temperature p|, |18||. Demanding the 



absence of such dangerous minima, the one loop analysis leads to an upper bound on 



phase transition strength. 
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the lightest CP-even Higgs mass of order 80 GeV. This bound was obtained for values 
of mfj = —mfj of order (80 GeV) 2 . 

The most important two loop corrections are of the form (ft 2 log(0) and, as said 
above, are induced by the Standard Model weak gauge bosons as well as by the stop 
and gluon loops P, 2"2]|. It was recently noticed that the coefficient of these terms can 



be efficiently obtained by the study of the three dimensional running mass of the scalar 
top and Higgs fields in the dimensionally reduced theory at high temperatures |30[ . 
Equivalently, in a four dimensional computation of the MSSM Higgs effective potential 
with a heavy left-handed stop, we obtain 
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where the first term comes from the Standard Model gauge boson-loop contributions, 
while the second and third terms come from the light supersymmetric particle loop 
contributions. The scale Ah depends on the finite corrections, which may be obtained 
by the expressions given in appendix A 0. As mentioned above, the two-loop corrections 
are very important and, as has been shown in Ref. [22], they can make the phase 
transition strongly first order even for mjj ~ 0. 

An analogous situation occurs in the [/-direction. The one-loop expression is ap- 
proximately given by 



V {U) + V 1 {U, T) = (-ml + 7c/ T 2 ) U 
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The contribution to Eu inside the squared brackets comes from the transverse gluons, 
while the one inside the curly brackets comes from the squark and Higgs contribu- 
tions |p~8 |. 

Analogous to the case of the field 0, the two loop corrections to the [/-potential are 
dominated by gluon and stop loops and are approximately given by 



V 2 (U,T) 
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where, as in the Higgs case, the scale Au may only be obtained after the finite corrections 
to the effective potential, given in appendix B, are computed. 



2 In the numerical computations, we use the whole expression given in appendix A. 
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Once the effective potential in the and U directions are computed, one can study 
the strength of the electroweak phase transition, as well as the presence of potential 
color breaking minima (see section 3). At one- loop, it was observed that requiring the 
stability of the physical vacuum at zero temperature was enough to assure the absolute 



stability of the potential at finite temperature. As has been first noticed in Ref. [30 
once two loop corrections are included, the situation is more complicated. We shall 
proceed with this analysis in the next section. 



3 Color breaking minima 

As explained above, we work in the limit of large values of uiq, for which small contribu- 



tions to Ap are expected. As has been shown in Ref. fllSfl , in this limit a color breaking 
minimum with (Q) ^ does not develop at T = 0, unless the stop mixing parameter 
A t is very large. These large values of A t suppress the strength of the first order phase 
transition and are hence of no interest for this study. We expect no modification of this 
conclusion at T ^ 0, so far as the relation ttiq 3> T is preserved. 

At zero temperature, for Q = the minimization of the effective potential for the 
fields and U shows that the true minima are located for vanishing values of one 
of the two fields. The two set of minima are connected through a family of saddle 
points for which both fields acquire non- vanishing values. Due to the nature of the high 
temperature corrections, we do not expect a modification of this conclusion at finite 
temperature. A detailed analysis of this problem would however be needed to decide 
whether minima with ^ and U ^ for Q = exist at finite temperature. We 
shall restrict our analysis to the behaviour of the electroweak and color breaking phase 
transitions in the directions ^ 0, U = and U ^ and = 0, respectively. 

Two parameters control the presence of color breaking minima: rnfj, defined as the 
smallest value of mjj for which a color breaking minimum deeper than the electroweak 
breaking minimum is present at T = 0, and , the critical temperature for the tran- 
sition into a color breaking minimum in the [/-direction. The value of may be 
obtained by analysing the effective potential for the field U at zero temperature, and it 
is approximately given by [TH|| 

m^p^)" 4 . (3.12) 

Defining the critical temperature as that at which the potential at the symmetry 
preserving and broken minima are degenerate, four situations can happen in the com- 
parison of the critical temperatures along the (T c ) and U (T^) transitions: 

a) T% < T c ; m v < m c v 



b) T? < T c ; mu > m L L 



c) > T c ; mu < m c v 

d) T% > T c ; mu > m c v 
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In the universe cools down, a phase transition into a color preserving 

minimum occurs, which remains stable until T = 0. This situation, of absolute stabil- 
ity of the physical vacuum, is the most conservative requirement to obtain electroweak 
baryogenesis. We shall generally demand these conditions to define the allowed pa- 
rameter space in the rrih-mj plane. In case b), at T = the color breaking minimum 
is deeper than the physical one implying that the color preserving minimum becomes 
unstable for finite values of the temperature, with T < T c . A physically acceptable 
situation may only occur if the lifetime of the physical vacuum is smaller than the age 
of the universe. In the following, we shall denote this situation as "metastability" . In 
case c), as the universe cools down, a color breaking minimum develops which, how- 
ever, becomes metastable as the temperature approaches zero. A physically acceptable 
situation can only take place if a two step phase transition occurs, that is if the color 
breaking minimum has a lifetime lower than the age of the universe at some temperature 



T < T c [ 50(1 . Finally, in case d) the color breaking minimum is absolutely stable and 
hence, the situation becomes physically unacceptable. 

As said before, at one-loop, the condition of having a color preserving minimum at 
zero temperature, mjj < mfj, automatically ensures that < T c . Hence, case c) never 
occurs and a two step phase transition is not allowed. On the contrary, as we shall show 
below, when two-loop corrections to the effective potential of the Higgs and stop fields 
are included, cases b) or c) may occur depending mainly on the values of the mixing in 
the stop sector and on tan (3. 

Figure 1 shows the region of parameter space consistent with a sufficiently strong 
phase transition for four different values of the stop mixing A t . At the left of the solid 
line, v/T > 1. Hence, for a given value of A t and rrih, the solid line gives the upper bound 
on the stop mass consistent with the preservation of the generated baryon number at the 
electroweak phase transition. Since v/T is inversely proportional to m\, lower values 
of mj (larger values of E) are needed as the Higgs mass is increased. The dashed line 
represents the values for which the critical temperature of the U field is equal to the 
critical temperature for the electroweak phase transition. This means that T c > T^ 
(T c < T^) on the right (left) of the dashed lines. Similarly, the short-dashed line is 
defined by the values of the Higgs and stop mass parameters for which my = m^, i.e. 
my < mfj (jhu > mfj) on the right (left) of the short-dashed lines. Observe that, for 
the same value of the stop mass parameters, larger Higgs masses are associated with 
larger values of tan (3. 

From Fig. 1, it is clear that for low values of the mixing, A t ;$ 200 GeV, case a) 
or c) may occur but, contrary to what happens at one-loop, case b) is not realized. 
For the case of no mixing, this result is in agreement with the analysis of [30]. The 
region of absolute stability of the physical vacuum for A t ~ is bounded to values 
of the Higgs mass of order 95 GeV. There is a small region at the right of the solid 
line, in which a two-step phase transition may take place, for values of the parameters 
which would lead to v/T < 1 for T = T c , but may evolve to larger values at some 
T < T c at which the second of the two step phase transition into the physical vacuum 
takes place. This region disappears for larger values of the stop mixing mass parameter. 
Quite generally, depending on the value of A t , the allowed region can be larger if a two 
step phase transition is considered, although the upper bound on the Higgs mass is not 
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modified (see section 4). For values of the mixing parameter A t between 200 GeV and 
300 GeV, both situations, cases b) and c) may occur, depending on the value of tan (3. 
For large values of the stop mixing, A t > 300 GeV, a two-step phase transition does not 
take place. The region of parameter space located at the left of both the dashed and 
short-dashed lines in the figure, corresponds to the case d) and hence is ruled out. 



4 Higgs and Stop mass constraints 

To define the allowed region of space in the mh-mj plane, one has to join all the allowed 
regions for different values of A t . As the value of A t increases, the stop mass values 
tend to decrease, and hence the lower bound on the stop mass is usually obtained 
for larger values of A t . Since for larger values of A t , the phase transition becomes 
weaker, large values of the Higgs mass may only be obtained for low values of A t . If 
we demand absolute stability, the upper bound on decreases quite fast with A t , and 
the Higgs mass is constrained to be below 105, 95, 85 GeV for A t = 300, 500, 600 GeV, 
respectively. Observe that, as becomes clear from Fig. 1, the same point in the rrih-mj 
plane may correspond to different values of the Higgs and stop parameters and hence 
could lead to different physical situations. 

Figure 2 shows the allowed region of parameter space in the m^-m^- plane. The 
solid line demarks the region where v/T > 1 and the physical vacuum is absolutely 
stable. The right handed boundary is defined by the points of v/T = 1 and A t ~ 
(values of A t > 100 GeV are necessary to reach the largest values of m h , see Fig. 1). 
The tiny region of parameter space, shown in the upper left panel of Fig. 1, for which 
a two step phase transition may occur for values of the stop mass at the right of the 
solid line, v(T c )/T c < 1, is not shown in the figure. Observe that the uppermost values 
of mh are approximately constant, independent of mr for a large region of parameter 
space. What constrains the Higgs mass in the region of A t ^ 300 GeV is not the 
mechanism of electroweak baryogenesis, but just the fact that, for the given values of 
A t and ttiq = 1 TeV, the Higgs mass cannot reach a larger bound. Slightly stronger 
(weaker) constraints would be obtained for rriq = 500 GeV (2 TeV). Finally, the left 
boundary is defined by values of rhu = mfj, for A t > 300 GeV. As said before, for this 
range of values of the stop mixing, and as happens at the one-loop level, the condition of 
stability at zero temperature ensures the stability at finite temperature, as it becomes 
clear by observing the evolution of the different regions in Fig. 1. 

All together, and even demanding absolute stability of the physical vacuum, elec- 
troweak baryogenesis seems to work for a wide region of Higgs and stop mass values. 



Higgs masses between the present experimental limit, of about 75 GeV pql , and around 
105 GeV are consistent with this scenario. Similarly, the running stop mass may vary 
from values of order 165 GeV (of the same order as the top quark mass one) and 100 
GeV. Observe that, due to the influence of the D-terms, values of 165 GeV, A t ~ 
and rrih ~ 75 GeV, are achieved for small positive values of ray (of order 40 GeV for 
the case represented in Fig. 2). Also observe that for lower values of itlq the phase 
transition may become more strongly first order and slightly larger values of the stop 
masses may be obtained. 
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The region bounded by the thin dashed lines in Fig. 2 corresponds to values of - 
m^-for which a two step phase transition can take place. Even though the effect of stop 
mixing greatly enlarges this region, we obtain from the present study that a two step 
phase transition may occur only in a narrow region of stop masses. The exact location 
of this band depends slightly on the value of tuq. For lower values of tuq, it moves to 
larger values of mj. Finally, the allowed parameter space may be greatly increased if 
we consider the condition of metastability. However, to analyse whether this may lead 
to a physically acceptable situation, the lifetime of the finite temperature metastable 
vacuum should be computed, something which is beyond the scope of this article. 



5 Conclusions 

In this article, we have computed the region of stop and Higgs masses consistent with the 
preservation of the baryon number generated at the electroweak phase transition. We 
showed that, once two loop corrections are consistently included in both the Higgs and 
the stop effective potentials, the electroweak phase transition may occur at one or two 
steps, with the additional possibility of a metastable physical vacuum state. Demanding 
the condition of absolute stability of the physical vacuum at all temperatures, we find 
that the phase transition may be sufficiently strong if the Higgs mass 75 GeV < < 
105 GeV and the lightest stop mass 100 GeV ^ mj ^ m t . The allowed m^-rrih region 
may be enlarged if the condition of metastability of the physical vacuum is considered, 
with stop mass parameters which are somewhat lower than in the case of an absolutely 
stable vacuum. A two step phase transition may only occur for values of the Higgs and 
stop masses which are also consistent with absolute stability, although the necessary 
stop mass parameters are slightly different. 

Our results present a strong case for Higgs and stop searches at LEP and the Teva- 
tron colliders. Indeed, by the end of 1999, LEP should be able to explore Higgs masses 
as high as 100 GeV, testing almost the whole region consistent with electroweak baryo- 
genesis. An upgrade in energy, up to \/s ~ 200 GeV during the years 1999-2000 will 
allow to explore the remaining region. If the Higgs is found at LEP, the Tevatron shall 
be able to test a wide region of the allowed stop masses during the next run, to start 
at the end of 1999. The stop mass reach of the upgraded Tevatron depends on the stop 
decay channels, but the whole region of stop masses consistent with this scenario will be 
explored if the ultimate upgraded Tevatron (TeV33) becomes operative. Observe that, 
LEP is beginning to test the lowest values of the stop mass |36| , which are compatible 
with a metastable electroweak vacuum state. 

Limits on the stop and Higgs mass parameters can also be obtained from elec- 
troweak precision measurements and rare decay processes. As explained in section 2, 
the parameter Ap sets a lower bound on the left handed stop mass. Another impor- 
tant measurement is the branching ratio BR(6 — > 57). Indeed, it has been recently 



argued [24] that the baryon number generation at the electroweak phase transition is 



suppressed for large values of the CP-odd Higgs mass []. Lower values of the CP-odd 

3 Observe, however, that both the Higgs mass and the value of v(T c )/T c vary only slightly with rriA, 
for values of tua ~ 250 GeV. The results presented in this article are hence expected to remain valid 
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Higgs mass, and hence of the charged Higgs mass, tend to enhance the branching ratio 
of this rare decay with respect to the Standard Model value |[7|. Since the Standard 



Model prediction for BR(6 — > 57) 0] is more than one standard deviation above the 
experimental value a consistent result may only be obtained if the stop-chargino 
contributions partially cancel the charged Higgs induced enhancement. For the val- 
ues of the parameters presented in Figure 1, and for a relatively light charged Higgs, 
rriH+ — 300 GeV (and hence ~ 300 GeV), this can only be achieved for moderate 
values of the stop mixing parameter A t (see ref. ||29|| ). The precise bound depends on 
the value of tan (3, /1 and the charged Higgs mass, but values of A t < 200 GeV tend to 
be disfavored. Although this fact does not modify the upper bound on the Higgs mass, 
it shows that the largest values of the stop masses, obtained for low values of A t tend 
to be disfavored (see Fig. 1). 

Finally, it is important to stress that, since the above results are obtained by using 
a two-loop finite temperature effective potential method, a non-perturbative analysis 
will be necessary to support these calculations, providing a definite check of the region 
of parameter space consistent with electroweak baryogenesis. 
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Appendix A 

In this appendix we present some technical details concerning the effective potential 
at finite temperature along the ^-direction, including the leading two-loop corrections. 
We will consider the case where the pseudoscalar mass is much larger than the 
temperature T. In this case only the light CP-even Higgs is coupled to the thermal bath, 
while the thermal distribution functions of the heavy neutral CP-even, the neutral CP- 
odd and the charged Higgses are Boltzmann suppressed and decouple from the thermal 
bath. We will work in the 't Hooft-Landau gauge and in the MS-renormalization scheme. 
We will follow the notation and conventions of Ref. |35] . Dimensional regularization 



(D = 4 — 2e) is used to evaluate divergent integrals. Poles in 1/e and terms depending 
on L e will be dropped since they are canceled by counterterms We will fix the 

MS-scale ~p to the temperature T, and hence put to zero all terms in ln(/l/T), and all 
the couplings and fields will be considered at the scale /I = T: g(T), 4>(T), ■ ■ ■. 

The states that give the most important contributions to the effective potential are 
the SM-particles: electroweak gauge bosons, the top quark, the Higgs and Goldstone 
bosons (W, Z, 7, t, h, x, respectively), as well as the left and right handed (third 
generation) scalar tops ( t^, £#). The SM fields have tree-level masses: 

2 1 2 1 2 



for these values of vha- 
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ml = \(g 2 + g' 2 ) <P 2 

m\ = ^ sin /3 2 h 2 <p 2 

ml = 1 -(g 2 + g> 2 )cos 2 2f3^ 2 -v 2 ) 

m 2 x = ±tf + g ») cos 2 2(3 (<p 2 - v 2 ) (A.l) 

and degrees of freedom (longitudinal, L, and transversal, T, for the gauge bosons) 

h Wl = 2 n Wr = 4 
nz L = 1 riz T = 2 

n h = 1 n x = 3 

= -12 . (A.2) 

The stop fields (ti, £r), with degrees of freedom, 

n T L = n T R = Q ( A - 3 ) 
have a squared mass matrix given by: 



M 2^ = I m 2 Q + m 2 t + \(g 2 - g 12 ) cos2p 4> 2 m t A t * 

' 1 m t A 4 mfj + m 2 + y 2 cos 2(3 (p 2 ' 



t \ <m . A . tt>2 — I — nm^ 1 _1_ —n f 2 



where the L — R mixing is defined as 

A t = A t -///tan/3 . (A. 5) 

The matrix ( |A.4| ) is diagonalized by the angle 6t defined by 



sin 20, = , 2M - cos 20, = , ^"-^ (A.6) 



(TrAl 2 ) 2 -4detA4 2 ^(TrAl 2 ) 2 - 4det A4 

and we will call t and T the light and heavy eigenstates, and mj and m ^ the corre- 
sponding eigenvalues. 

We now define thermal masses for the longitudinal components of gauge bosons 
(Wl, Zl,Jl) and for all scalars (h,x, tL, £r) by means of the corresponding self- 
energies. We include in the thermal bath the SM particles (W, Z, 7, t, h and x)> 
as well as the light stop t, which plays a prominent role in the strength of the phase 
transition, and (light) charginos and neutralinos, which were proven to be essential for 



the generation of the baryon asymmetry p4| . The self-energies, to leading order, are 
given by: 

1% = \g 2 T 2 
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n 



B 



-g' 2 T 2 
9 y 



U h = — (g 2 + g' 2 ) cos 2 2(3 T 2 + —g 2 T 2 + —g' 2 T 2 + -h 2 sin 2 (3 T 2 
16 v ; 16 48 2 * 



n 



n 



n 



^ T 2 + T 2 - -Lg' 2 (1 + cos 2/3) T 2 + \-h\ (2 + sin 2 /?) T 2 



1^2 T 2 1 ( g 2/3) ,2 T 2 1^2 

9 s 18 v 7 6 ' 



l + sm 2 /3(l-4 

\ TTln 



T 2 . (A.7) 



The thermal masses are defined, for the SM bosons, as: 



™2 



??? 



1 

2 



4 



(g' + g'^ + Uw + IlB 



± 



\ 



m 2 



and for the stops as: 



(A.f 



m g + m t + 1(9 - 9 ) cos 2/5 2 + n Ti 



m t A t m 2 / + m 2 + |g' 2 cos 2/? (b 2 + II ^ 



(A.9) 

Similarly one can diagonalize the mass matrix (|A.9|) by means of the angle 9 t , leading 
to the mass eigenstates m^-and m ^. For mq 3> T: 



sin 6+ ~ sin & 



(A.10) 



We can now expand the effective potential as a sum 

V(<b,T) = V + V 1 + V 2 + -- 



(A.ll) 



where V n indicates the n-th loop potential in the resummed theory at finite temperature. 
The tree-level potential can be written as: 



V = -\m 2 ^)cb 2 + ^(9 2 + 9' 2 ) cos 2 2(3 4 



(A.12) 



where 



j ji ■ 

m 2 (jx) = -m 2 z (v) cos 2 2(3 + J2 J^^^ 



dmfiv) 



dv 2 



log(m 2 (v)/fi 2 ) + 7;~Ci 



(A13) 
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where i = W, Z, h, x, t, t, T, the constants are Cj = 5/6 (Cj = 3/2) for gauge bosons 
(scalar bosons and fermions), and the VEV of the Higgs field is normalized to v = 246 
GeV. 

We will perform daisy resummation on the n = modes of the longitudinal compo- 
nents of the gauge bosons Wl, Zl, 7l and for all modes of the bosons h, x and the light 
stop t. No resummation is performed on fermions (i.e. the top quark) or for heavy 
bosons which decouple from the thermal bath (in particular the heavy stop T). With 
the above prescription we can write the one-loop effective potential V\ as: 



64tt 2 



(A.14) 



where, as said before, for the analysis of the phase transition, we shall take p, = T. In 
the above, the sum extends over i = Wl, Zl, Jl, Wt, Zt, t, h, x, t, T, and the masses M; 
are defined by 

rrii i = W L ,Z L ,j L ,W T ,Z T ,t,f 



(A.15) 



mi i = h, x, t . 
The thermal contributions J® are defined by 



j(0 



Mm!) 
Jr(tT&) 



mj — m- 



t J F (m 2 i) 
and the thermal integrals Jb,f are: 

JbAv 2 ) - 



W l ,Z l ,1l 

W T ,Z T ,f 

h,xJ 
t , 



(A.16) 



dx x 2 log If e 



(A.17) 



The two-loop potential can be split into the SM (V 2 SM ) and the MSSM (\/ 2 MSSM ) 
terms. 

V 2 = V 2 SM + K 2 MSSM (A. 18) 

The value of V 2 SM has been computed in Ref. |35|] where the full expression can be found. 
It is dominated by logarithmic terms that can be written as: 



SM 



9 _rp2 



16tt 2 



3 -m 2 w + 2U W log 



m w + 2m WL 63 
3T ~8 



2 , m W 



3 2 , 2m Wr 
+ -m™ log — 

4 3T 



3 2 2m w Q _ 



(A.19) 



where the approximation g' = has been used. There are also non-logarithmic contri- 
butions involving the strong g s and Yukawa h t couplings which are less important for 
the phase transition. 
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In the two-loop potential y^ MSSM the strong and top Yukawa couplings are involved 
in logarithmic contributions. Therefore we have found to be an excellent approximation 
to keep g s and ht, and to neglect the weak couplings g and g'. In this approximation the 
relevant diagrams that contribute to V^ MSSM are of two kinds: sunset diagrams, labeled 
by Vxyz, where X, Y and Z are the propagating fields, and figure eight diagrams, 
labeled by Vxy, with propagating X and Y fields. With this prescription can 
be decomposed as: 

yMSSM = y_ + y m + y^ + y^ + y^ + y_ ^ ^ 

where g stands for gluons. The different contributions are given by: 

V~ = - g l(Nl-i)v ssv {m Tl m T ^) 

Vjj h = — (h 2 sin 2 p -^j= + h t A t sin f3 sin 9 t cos OA N c H(m h) m^, m~) T 2 



9 

1 



2 



V g j = -^{Nt-l)V sv {m Tl U) 



V n = ^h 2 t sm 2 f3N c I{m T )I{m h ) 
3 



V 7x = "-h 2 sm 2 pN c I(jn~)I(jn x ) 
V tt = 4 N c (N c + l)I 2 (mj) (A.21) 



The functions involved in ( [A..21 ) are all defined in Ref. [3~5|. We can use them with the 
following prescriptions: 

• Cancel all poles in 1/e, and i e -terms, with counterterms. 

• Take the limit e — > and keep all finite contributions. 

• Cancel all cb terms for scalar bosons with thermal counterterms. 

• Linear terms as MjT 3 which appear in figure-eight diagrams are canceled by ther- 
mal counterterms. 

Following these rules, we obtain results which are in good agreement with those 



found in Ref. [22 



Appendix B 

In this appendix we will specify the effective potential at finite temperature in the 
background field U = t^Ua, where u a is a constant unit vector in color space Q, which 
breaks SU C (3) into 577(2). We will proceed as in appendix A and present the result of 

4 Using the SU C (3) invariance the vector u a can be chosen as u a = (1,0, 0). 
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the two-loop calculation. We will neglect g' and keep the strong g s and top Yukawa h t 
couplings. In the following, we use g' to denote one of the gluon states appearing in the 
spectrum. 

The states contributing to the effective potential |1| are the four gluons g and the 
gluon g', five real squarks uj (would-be goldstones) and the real squark p, four light Ti 
and four heavy Q scalars, coming from the mixing between the SM Higgs doublet H and 
the left-handed (third generation) squark doublet q = q a u Q , and two degenerate massive 
Dirac fermions / coming from the mixing between the left-handed (third generation) 
fermion doublet q^ = q1u a and the higgsino fl8|, |3(J. The masses of g, g', u, p and / 
are given by: 



m 2 g , 



l& 2 



and their degrees of freedom are 





2 , 1 2 7-7-2 


m 2 


2 i 2rr2 


m) 


= p 2 + h 2 t U 2 


n 9L 


= 4 n 9T = 8 


U 9' L 


= 1 ^ = 2 




= 5 n p = 1 


n f 


= -8 . 



(B.1) 



The scalars (H, q), with degrees of freedom 

n H = 4 n~ = 4 

have a squared mass matrix 

' -\m\ + h\ sin 2 /? f/ 2 h t sin /3 f/l t 



^,7 



ht sin /3 C/A t 



m| + h 2 U 2 



(B.2) 



(B.3) 



(B.4) 



After diagonalization of the matrix ( p.4j ) by the rotation defined through the angle 6-n, 
defined as in ( |A.6| ), we obtain the light Tt and heavy Q eigenstates and the corresponding 
eigenvalues, m 2 ^ and m 2 Q . 

Using the self-energies (|A.7|) , and the gluon self-energy 



-g 2 T 2 
3 ys 



(B.5) 



we can define the thermal masses for the longitudinal components of gauge bosons (gi, 
g' L ) and for all scalar bosons (a;, p, H, Q) as: 



m 9L + n g 
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m 9' L 

m 2 



m g' L + U 9 



(B.6) 



while and mg are the eigenvalues of the thermally corrected squared mass matrix 

/ 1™2 , u2 ;„2 fl rr2 



Al 



£T, g 



2 -m£ + /i 2 sin 2 (3 U 2 + Il h h t sin (3 UA t 
h t sin (3 UA t m 2 Q + h 2 U 2 + U 



t L 



diagonalized by the rotation H . For mq ^> T: 



sin dfi ~ sin 



h t sin (3 UA t 



m. 



We will now expand the effective potential in a loop expansion as in (|A.ll ), 

V{U,T) = V + V 1 + V 2 + ... 
The tree level potential can be written as: 



V = mlifi) U 2 + -g 2 U 4 



with 



2 ,_s 2 Ui 2/ \dTn 2 (u) 



32tt 



dv? 



io g K 2 ( M )//i 2 ) + - - Ct 



(B.7) 

(B.8) 

(B.9) 
(B.fO) 

(B.H) 



where, as explained in Appendix A, we shall take \i = T . In the above, i = g, g', lu, p, f, H, Q 



and u = (U) = ^ -3m 2 ; / ' g 2 s . 

We will perform, as in the case of appendix A, resummation on the n = modes 
of the longitudinal components of the gauge bosons gL^g'i an d for all modes of the 
bosons u,p and the light scalar 7i. No resummation is performed either on fermions / 
or for the heavy scalar boson Q, whose mass is dominated by the soft-term itlq and is 
decoupled from the thermal bath. We can then write the one-loop effective potential 
as in Eq. ( |A14| ), where the sum extends over i = g L , g' L , g T .g' T , f,u, p,H, Q, and the 
masses M,- are defined as: 



rrii 



gL,g' L ,gT,g' T J, Q 



Ma 



;b.i2) 



rrii i = u, p,H 
The thermal contributions jW are defined 



JB(m 2 ) - - (mf - 

J B (m 2 ) 
J B {m 2 ) 
J F (m 2 ) 



rri; 



9l, 9l, 1l 
9t,9t, Q 

f ■ 



(B.13) 
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The two loop diagrams that contribute to V-i are, as in appendix A, of two kinds: 
sunset diagrams labeled by Vxyz, where X, Y and Z are propagating fields, and figure 
eight diagrams labeled by Vxy, with propagating X and Y fields. In the following 
we will denote: G = (g,g r ) and t = (uj,p), rj being the ghost fields. The relevant 
contributions to V2 are then (see also Ref. |3( 



V< 



g<;g = ~g 2 s ^ [(N c - 2)V vvv (m g ,m g ,0) + V vvv (m g ,m g ,m g/ )] 



VvvG = -9s^T \- 2 ( N c ~ l)^v(0, 0, m g ) + V vvV (0, 0, m g , 



V ttc = i( N c ~ ^Pssv{m u , m w , m g ) + (N c - l)£> 5S y(m w , m p , m g ) 

+ N ° N 1 ^ssv(jn ul ,m p ,m g/ ) + ^Vssv(?f^,rfi u ,™>g>) 
+ N e (N c -2)V ssv (m u ,m u ,0)] 



V 



tGG 



2 m l 



_ (jv c -l) 2 _ 
(JV C - l)V svv (m p , m g , m g ) + 2 — — V S vv{m p , m g >, m g >) 



+ 



(Nc - 2J 

N c 

,N C 



Vswirriu, m g , m g >) + N C (N C - 2)V svv (m U) , m g , 0) 



V GG = -9t-£ [2(A^ C - 2)V VV (0, m g ) + 2V vv (m g , m g ,) + (N c - l)V vv (m g , m g )\ 



V 



tG 



= --^{(Nc- 1 ) [3P 5 y(^, mg) + V sv (m p , m g )} 

o 

+ y ^ Nc + l ) V sv{Tn^ m g ,) + (N c - l)V sv (m p , m g ,)} 
+ 2N c (N c -2)V sv (m u ,,0)} 



V 



1 1 1 



V 



tnn 



V 



1 1 



-f§ [3H(m p , m p , m p ) + (2N C - l)H(m p , m u , mj T 2 U 2 

io 

—2 h 2 sin 2 /3 U + h t A t sin j3 sin 6*^ cos 9<n Hijn^^m^m-^T 2 
^g 2 [3I 2 (m p ) + (AN C - 2)I(jn p )I(jn u ) + (AN 2 - l)I 2 (m w ) 



V In = h 2 t sin 2 p cos 2 9 H I(m H ) [I(m p ) + (2N C - 1)/ 



(B.14) 



where all functions involved in ( B.14j ) are defined in Ref. and cancellation of poles 
in 1/e, and i e and cb terms, as well as linear terms in MjT 3 , is achieved as in the case 
of appendix A against counterterms and thermal counterterms. 
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Figure 1: Values of m h , mj for which v(T c )/T c = 1 (solid line), = T c (dashed line), 
rhu = rnfj (short-dashed line), for tuq = 1 TeV and different fixed values of A t . The 
region on the left of the solid line is consistent with a strongly first order phase transition. 
A two step phase transition may occur in the regions on the left of the dashed line, while 
on the left of the short-dashed line, the physical vacuum at T = becomes metastable. 
The region on the left of both the dashed and short-dashed lines leads to a stable color 
breaking vacuum state at zero temperature and is hence physically unacceptable. 
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mr (GeV) 

Figure 2: Region of the m/j-m j parameter space for which a strongly first order phase 
transition takes place is shown within solid lines. The short-dashed lines demark the 
region for which a two-step phase transition may occur. The region on the right of the 
dashed line and left of the short-dashed may lead to a metastable vacuum state. 
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